Spin current and rectification in Luttinger liquids 
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We demonstrate that spin current can be generated by an ac voltage in a one-channel quantum 
wire with strong repulsive electron interactions in the presence of a non-magnetic impurity and 
uniform static magnetic field. In a certain range of voltages, the spin current can exhibit a power 
dependence on the ac voltage bias with a negative exponent. The spin current expressed in units 
of ft/2 per second can become much larger than the charge current in units of the electron charge 
per second. The spin current generation requires neither spin-polarized particle injection nor time- 
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I. INTRODUCTION 

The pioneering paper by Christen and Biittikcr 1 has 
stimulated much interest to rectification in quantum 
wires and other mesoscopic systems. Most attention was 
focused on the simplest case of Fermi liquids 2 -^. Recently 
this research was extended to strongly interacting sys- 
tems where Luttinger liquids are formed4^S^. One of 
the topics of current interest is the rectification effect in 
Luttinger liquids in a magnetic fi.eld3 » &£. In the presence 
of a magnetic field, both spin and charge currents can be 
generated. So far, however, only charge currents in Lut- 
tinger liquids have been studied. In this paper we show 
that a dc spin current can be generated by an ac voltage 
bias in a single-channel quantum wire. 

In recent years many approaches to the generation of 
spin currents in quantum wires were put forward. Typ- 
ically both, spin and charge currents are generated and 
the spin current expressed in units of h/2 per second is 
smaller than the electric current in units of e per sec- 
ond (e is the electron charge). Such a situation naturally 
emerges in partially polarized systems since each electron 
carries the charge e and its spin projection on the z-axis 
is ±H/2. A proposal how to obtain a spin current exceed- 
ing the charge current in a quantum wire was published 
by Sharma and Chamon^ who considered a Luttinger liq- 
uid in the presence of a time-dependent magnetic field in 
a region of the size of an electron wavelength. In a very 
different physical context, a spin current without charge 
current was predicted for edge modes in the quantum 
Hall effect in graphene 9 . Pure spin currents can also flow 
in open circuits which cannot support charge currents— 
In this paper, we show that the generation of a dc spin 
current exceeding the charge current is also possible in 
closed circuits without time-dependent magnetic fields. 
The spin current can be generated in a spatially asym- 
metric Luttinger liquid system in the presence of an ac 
bias. Interestingly, in a certain interval of low voltages 
the dc spin current grows as a negative power of the ac 
voltage when the voltage decreases. 

The paper is organized as follows: The next section 



contains a qualitative discussion. We briefly address the 
simplest case of non-interacting electrons, discuss its dif- 
ferences from the most interesting case of strong elec- 
tron interaction and estimate at what conditions the ef- 
fect can be observed. Section IIIII contains the details 
of the bosonization procedure which we use to treat the 
electron-electron interaction. In Section [TVl we calculate 
analytically the spin and charge rectification currents in 
the presence of a weak asymmetric potential. Numerical 
results for a simple model with strong asymmetric poten- 
tial are discussed in Appendix [A] Appendices [5] and [C] 
contain technical details of the perturbation theory em- 
ployed in Sec. [IV] 



II. MODEL AND PHYSICS OF THE PROBLEM 

The rectifying quantum wire is sketched in Fig. [T] It 
consists of a one-dimensional conductor with a scatterer 
in the center of the system at x = 0. The scatterer 
creates an asymmetric potential U(x) ^ U{—x). The size 
of the scatterer ajj ~ l/kp is of the order of the electron 
wavelength. A spin current can be generated only if time- 
reversal symmetry is broken. Thus, we assume that the 
system is placed in a uniform magnetic field H. The 
field defines the S z direction of the electron spins. If the 
wire is sufficiently narrow then the effect of the magnetic 
field on the kinetic energy of electrons can be neglected 
and the field enters the problem only via its interaction 
with the spins. At its two ends, the wire is connected 
to nonmagnetic electrodes, labeled by i = 1, 2. The left 
electrode, i = 1, is controlled by an ac voltage source, 
while the right electrode, i = 2, is kept on ground. 

The magnetic field H breaks the symmetry between 
the two orientations of the electron spin. In a uniform 
wire this would not result in a net spin current since the 
conductances of the spin-up and -down channels would 
be the same^ e 2 /h, and the spin currents of the spin- 
up and -down electrons would be opposite. In the pres- 
ence of a potential barrier such a cancellation does not 
occur—. In a system with strong electron interaction, 
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the spatial asymmetry of the wire leads to an asymmet- 
ric I — V curved I(V) ^ I{—V). Thus, an ac voltage 
bias generates spin and charge dc currents, and I r c . 

As we will see, the problem is most interesting in the 
case of strong electron interaction. Before addressing 
that more difficult case, let us discuss what happens in 
the absence of electron interaction. By non-interacting 
system we mean a wire in which electron-electron inter- 
action is completely screened by the gates. In such sit- 
uation the charge density in the wire is not fixed but 
depends on the gate potential and the electrochemical 
potentials of the leads. The leads define the chemical 
potentials [ip and /ip of the left- and right-moving elec- 
trons which are injected from the right and left reser- 
voirs respectively. In what follows we will assume that 
the chemical potential /xl = and [Ir oscillates between 
+eV and —eV. Thus, in the presence of the magnetic 
field H, the Fermi energies counted from the band bot- 
tom equal Ep(S z ) = Ep + 2S z /j,H for the left-moving 
electrons and Ep(S z ) = Ep + fiR + 2S z fiH for the right- 
movers, where S z = ±1/2 is the electron spin projection, 
Ep the Fermi level in the absence of the magnetic field 
and voltage bias and /j, the electron magnetic moment. 

As we will see, in a strongly interacting system the 
form of the potential barrier U(x) plays little role. How- 
ever, it is crucial in the non-interacting case. Let us 
choose U(x) in the form of the double potential barrier 
so that quasistationary levels E n , n = 0, 1, . . . are present 
(Fig. [2|). Thus, in the non-interacting case we consider a 
resonant tunneling diode 11 . The spin and charge currents 
as the functions of the chemical potential /j,r are 

I c { m ) = I(fi R , 1/2) + 7(/i fl , -1/2); (1) 

him) = 1/2) - -i/2)], (2) 

Ze 

where e is the electron charge, 

I(»R,S Z ) = I j " * dET(E); (3) 

and T(E) is the transmission coefficient. The transmis- 
sion coefficient is small far from the energies of the quasis- 
tationary levels E = E n and increases as E approaches 
E n . Obviously, if one applies a dc voltage eV = lip 
then the charge current, expressed in units of e per sec- 
ond, exceeds the spin current in units of h/2 per second. 
The situation changes in the presence of an ac bias. The 
dc currents generated by an ac voltage bias can be es- 
timated as I r . a — [I c / S {y) + I c / s (—V)]/2. Let us now 
assume that the magnetic field is tuned such that the 
Fermi level of the spin-down electrons Ep — fiH is close 
to the quasistationary level Eq and exceeds Eq while the 
Fermi level of the spin-up electrons Ep + /j,H is close 
to Ei and lies below E\. Let also eV be smaller than 
the distances, \E — Ep + fiH\ and \Ex — Ep — [iH\, 
between the Fermi levels and quasistationary levels in 



the absence of the voltage bias. In addition, we assume 
that T(E F + fiH) = T{E F - fiH). From the energy de- 
pendence of the transmission coefficient T(E) near the 
resonant levels, one finds that 7(eV,l/2) > I(eV, — 1/2) 
and \I(-eV, 1/2) | < |I(-eV,-l/2)|. Hence, F s > ^I r c . 
By an appropriate choice of parameters, one can produce 
any ratio of the spin and charge rectification currents. 

Note that the rectification effect for non-interacting 
electrons is possible even if the potential U{x) is sym- 
metric. The asymmetry of the system, necessary for rec- 
tification, is introduced by the applied voltage bias. The 
charge density injected into the wire from the leads is 
proportional to /zj, + [ip and hence is different for the 
opposite signs of the voltage. If the injected charge den- 
sity were independent of the voltage sign, i.e. )jlr os- 
cillated between eV/2 and — eV/2 and [ip = — hr os- 
cillated between —eV/2 and eV/2, then the rectification 
effect would be impossible for non-interacting electrons. 
This follows from Eqs. (JTJSJ) and the fact that the trans- 
mission coefficient T(E) is independent of the direction 
of the incoming wave for non- interacting particles 12 . In 
the presence of electron repulsion both the asymmetry of 
the potential and the voltage dependence of the injected 
charge contribute to the rectification current. It turns 
out, that in the case of strong electron interaction, the 
rectification effect due to the asymmetry of the potential 
barrier dominates. 

The above example is based on a special form of the 
potential barrier in the wire and assumes that the mag- 
netic field and chemical potentials are tuned in order to 
obtain the desired effect. As shown below, in the pres- 
ence of strong repulsive electron interaction no tuning 
is necessary and no quasistationary states are needed to 
obtain the spin current which is greater than the charge 
current. In fact, the spin rectification effect is possible 
even for weak asymmetric potentials U{x). This can 
be understood from the following toy model (a related 
model for rectification in a two-dimensional electron gas 
was studied in Ref . HH) : Let there be no uniform mag- 
netic field H and no asymmetric potential U(x). Instead, 
both right<-s-left and spin-up^->spin-down symmetries are 
broken by a weak coordinate-dependent magnetic field 
B z (x) ^ B z (—x), which is localized in a small region of 
size ~ 1/kp (we do not include the components B x y in 
the toy model) . Let us also assume that the spin-up and 
-down electrons do not interact with the electrons of the 
opposite spin. Then the system can be described as the 
combination of two spin-polarized one-channel wires with 
opposite spin-dependent potentials ±/j,B z (x), where /z is 
the electron magnetic moment. According to Ref. an ac 
bias generates a rectification current in each of those two 
systems and the currents are proportional to the cubes of 
the potentials (±/j,B z ) 3 . Thus, = — JJ. Hence, no net 
charge current I r = 7J + 7J is generated in the leading 
order. At the same time, there is a nonzero spin current 
in the third order in B z . A similar effect is present in a 
more realistic Luttinger liquid model considered below. 

The main focus of this paper is on the case of weak 
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asymmetric potentials. A simple model with strong im- 
purities is studied in Appendix [A] In Figs. [3] and H] we 
have represented the results from a numerical evaluation 
for the spin and charge currents I r s c for the potential 
shown in Fig. [2j Fig. [3] shows the non-interacting case. 
In Fig. |4]we represent the case of strong electron interac- 
tion. We have chosen parameters (explained in the figure 
captions) such that l r c is smaller than JJ for a range of 
the applied voltage. Further information on the numeri- 
cal approach is given in Appendix lAl 

Transport in a strongly interacting system in the pres- 
ence of a strong asymmetric potential U{x) is a diffi- 
cult problem which cannot be solved analytically and 
is sensitive to a particular choice of the potential. As 
we have mentioned, Appendix [Al contains the numerical 
analysis of a simple model of interacting electrons with 
a strong potential barrier. On the other hand, the inter- 
acting problem can be solved analytically in the limit of 
a weak potential U (x) with the help of the bosonization 
and Keldysh techniques (Sec. IIII[) . We will see that the 
rectification current exhibits a number of universal fea- 
tures, independent of the form of the potential U(x). In 
particular, in a wide interval of interaction strength, the 
spin rectification current can exceed the charge rectifi- 
cation current for an arbitrary shape of the asymmetric 
potential barrier. 

Rectification is a nonlinear transport phenomenon. 
Thus, it cannot be observed at low voltages at which the 
I — V curve is linear and hence symmetric. In Luttinger 
liquids the I— V curve is nonlinear at eV > ksT, where T 
is the temperature^. We will concentrate on the limit of 
the zero temperature which corresponds to the strongest 
rectification. We expect qualitatively the same behavior 
at T ~ V. At higher temperatures the charge and spin 
rectification effects disappear. Since the temperatures of 
the order of millikelvins can be achieved with dilution 
refrigeration, the rectification effect is possible even for 
the voltages as low as V < 1 /zV. 

In this paper we focus on the low- frequency ac bias. 
We define the rectification current as the dc response to 
a low-frequency square voltage wave of amplitude V: 

r s (V) = [I s (V)+I s (-V)}/2, (4) 

F c (V) = [I c (V)+I c (-V)}/2. (5) 

The above dc currents express via the currents of spin-up 
and -down electrons: F c = J[ + ij, F s = (h/2e)[I\ - /[]. 
The spin current exceeds the charge current if the signs 
of ij and 7[ are opposite. Equations (|4I5[) for the dc- 
currents do not contain the frequency uj of the ac-bias. 
They are valid as long as the frequency 

w < eV/h. (6) 

Indeed, as shown below, the rectification current is de- 
termined by electron backscattering off the asymmetric 
potential. Hence, one can neglect the time-dependence 
of the ac voltage in Eqs. (|4l5p if the period of the ac 



bias exceeds the duration r of one backscattering event. 
The time r ~ Tt ravo i + T uncer t a inty includes two contribu- 
tions. Travel is the time of the electron travel across the 
potential barrier. T uncerta inty comes from the uncertainty 
of the energy of the backscattered particle. If the barrier 
amplitude U(x) < Ep and the barrier occupies a region 
of size a J/ ~ 1/kp then Tt ravc i ~ 1/kpv ~ h/Ep, where 
v ~ hkp/m is the electron velocity. The energy uncer- 
tainty ~ eV translates into r unC ertainty ~ h/eV. Thus, 
for eV < Ep one obtains the condition ©. The same 
condition can be derived with the approach of Appendix 
A of Ref. [l5| and emerges in a related problem^. Note 
that for realistic voltages the low-frequency condition ([6]) 
allows rather high frequencies. Even for V ~ 1 fiY the 
maximal u) *~ 1 GHz. 

There remains the question of the asymmetric im- 
purity: We require a potential U(x) that is localized 
within ~ 1/kp. A possible realization is to generate 
two different (symmetric) local potentials by two gates 
within a distance ~ 1/kp or an electric potential cre- 
ated by an asymmetric gate of size ~ 1/kp placed at 
the distance ~ 1/kp from the wire. Electron densities 
of p ~ 10 11 cm -2 are possible nowadays in 2-dimcnsional 
electron gases, yielding 1/kp up to several 10 nm. Con- 
finement in a one-dimensional wire will reduce the elec- 
tron density further so that this number may increase 
further. Modern techniques allow placing electric gates 
of widths of ~ 20 nm at distances of ~ 20 — 50 nm. 
A realization of an asymmetric potential in this way is, 
therefore, within the reach. Alternatively, in the case 
of shorter electron wave-length, it should be possible to 
place an asymmetrically shaped STM tip close to the 
wire. An applied bias would yield an asymmetric scat- 
tering potential. With such a tip the asymmetry cannot 
be directly tuned, but most of our predictions are not 
sensitive to the precise shape of the potential. Certainly, 
an asymmetric potential may simply emerge by chance 
due to the presence of two point impurities of unequal 
strength at the distance ~ 1/kp. 



III. BOSONIZATION AND KELDYSH 
TECHNIQUE 

At w < V, the calculation of the rectification cur- 
rents reduces to the calculation of the stationary con- 
tributions to the dc I — V curves I S (V) and I C (V) that 
are even in the voltage V. We assume that the Coulomb 
interaction between distant charges is screened by the 
gates. This will allow us to use the standard Tomonaga- 
Luttinger model with short range interactions^ Elec- 
tric fields of external charges are also assumed to be 
screened. Thus, the applied voltage reveals itself only 
as the difference of the electrochemical potentials E\ and 
£2 °f the particles injected from the left and right reser- 
voirs. We assume that one lead is connected to the 
ground so that its electrochemical potential E2 = Ep 
is fixed. The electrochemical potential of the second lead 
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Ei = Ep + eV is controlled by the voltage source (see 
Fig. [TJ. Since the Tomonaga-Luttinger model captures 
only low-energy physics, we assume that eV < Ep, where 
Ep is of the order of the bandwidth. Rectification occurs 
due to backscattering off the asymmetric potential U(x). 
We will assume that the asymmetric potential is weak, 
U (x) < Ep. This will enable us to use perturbation the- 
ory. 

We assume that the magnetic field H couples only to 
the electron spin and we neglect the correction — eA/c 
to the momentum in the electron kinetic energy. Indeed, 
for a uniform field one can choose A ~ y, where the 
y-axis is orthogonal to the wire, and y is small inside a 
narrow wire. As shown in Ref.[l8|, such a system allows a 
formulation within the bosonization language and, in the 
absence of the asymmetric potential, can be described by 
a quadratic bosonic Hamiltonian 



*o= E E 

v.y' — L,R a.a' — ] 



dx{3 x (j) V! y)T-L V aya' {d x 4>y'a'), 



where a is the spin projection and v = R, L labels the 
left and right moving electrons, which are related to the 
boson fields <f> va as ipt a (x) ~ ^ v<je ±<k Fua x+ct> ua {x)) 
i for v = L. The operators are the Klein factors 
adding a particle of type {v, a) to the system, and kp^^/ir 
is the density of (v, a) particles in the system. The den- 
sities of the spin-up and -down electrons are different 
since the system is polarized by the external magnetic 
field. The 4x4 matrix Ti describes the electron-electron 
interactions. In the absence of spin-orbit interactions, 
L <-> R parity is conserved and we can introduce the 
quantities 4> a = 4> La + 4>R a and II CT = 4> La - 4> Ra such 
that the Hamiltonian decouples into two terms depend- 
ing on <f> a and only. In the absence of the external 
field, this Hamiltonian would further be diagonalized by 
the combinations 4> CjS oc </>j ± (j>i, and similarly for n c jS , 
expressing the spin and charge separation. This is here 
no longer the case because of the external magnetic field. 
If we focus on the <fi fields only (as n will not appear 
in the operators describing backscattering off U(x)), the 

, have a more 
which we can write 



fields diagonalizing the Hamiltonian, 



complicated linear relation to 
as 



a] 



-0] 



(8) 



and which corresponds to the matrix A T of Ref. 0. The 
normalization has been chosen such that the propagator 
of the <f> fields with respect to the Hamiltonian jTJ evalu- 
ates to {4> c ,s(ti)4>c,s(t2)) = -21n(i(£i ~t 2 )/T c + 5), where 
8 > is an infinitesimal quantity and t c ~ h/ Ep the ul- 
traviolet cutoff time. For non-interacting electrons with- 
out a magnetic field, g c — g s = 1/2. g c < 1/2 (> 1/2) for 
repulsive (attractive) interactions. The interaction con- 
stants depend on microscopic details and the magnetic 
field. The dimensionless parameter which controls the 



interaction strength is the ratio of the potential and ki- 
netic energies of the electrons. This ratio grows as the 
charge density decreases and hence lower electron densi- 
ties correspond to stronger repulsive interaction. In the 
absence of the magnetic field, terms in ([7]) in the form 
of exp(±2iy^70 s ) may become relevant and open a spin 
gap for g s < 1/2. In our model they can be neglected 
since they are suppressed by the rapidly oscillating fac- 
tors exp(±2i[fcj?f — kpi\x). It is convenient to model the 
leads as the regions near the right and left ends of the 
wire without electron interaction^. 

Backscattering off the impurity potential U(x) 
is described by the following contribution to the 
Hamiltonianii H = Hn + H': 



H'= E U(n hni )e 



m T T (O)+irii0|(O) 



(9) 



where the fields are evaluated at the impurity position 
x = and U(n-\, ni) = U*{— rtf, —ni) since the Hamilto- 
nian is Hermitian. The fields n do not enter the above 
equation due to the conservation of the electric charge 
and the z-projection of the spin. The Klein factors are 
not written because they drop out in the perturbative 
expansion. U (n j , n j ) are the amplitudes of backscat- 
tering of n-f spin-up and spin-down particles with 
n a > for L — > R and n a < for R — > L scat- 
tering. U(n^,ni) can be estimated as^ U(n-\,ni) ~ 
k F J dxU{x)e m i 2kF i x+m t 2kF ± x - U, where U is the max- 
imum of U{x). In the case of a symmetric potential, 
U(x) = U(—x), the coefficients £/(rif,n|) are real. 

The spin and charge current can be expressed as 



I s ,c — lI c 



nl _ r2 



Rr. 



(10) 



where L\ c and R\ c denote the current of the left- and 
right-movers near electrode i, respectively (see Fig. 
For a clean system (U(x) = 0), the currents obey 1 ? 
Rl c = fl2 c) L i c = L 2 c and Ic = 2e 2 V/h, I s = 0. With 

backscattering off U (x), particles are transferred between 
L and R in the wire, and hence R 2 S = R\ + dSn/dt, 
R 2 ': = R\ + dQu/dt, where Qr and Sr denote the to- 
tal charge and the ^-projection of the spin of the right- 
moving electrons^. The currents L 2 c and R} s c are de- 
termined by the leads (i.e. the regions without electron 
interaction in our modeli^) and remain the same as in the 
absence of the asymmetric potential. Thus, the spin and 
charge current can be represented as I c = 2e 2 V/h + 
and I s = I* s , where the backscattering current operators 
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i b c s =dQ R /dt = i[H,Q R ]/h 

= IT S ("T + ^)C/(«T^l)e mT0T(O)+ ™ i0i(O) , 

(ii) 

7 s 6s = d&i/dfc 

= -| E ("T -n i )C/(n T) n i )e i ™t^T(o)+i«i^(o). 

(12) 

The calculation of the rectification currents reduces to 
the calculation of the currents (fTTj) . (fT2|) at two opposite 
values of the dc voltage. 

To find the backscattered current we use the Keldysh 
technique^. We assume that at t — — oo there is no 
backscattcring in the Hamiltonian (U(x) — 0), and then 
the backscattering is gradually turned on. Thus, at 
t = — oo, the numbers Nl and Nr of the left- and right- 
moving electrons conserve separately: The system can 
be described by a partition function with two chemical 
potentials E\ = Ep + eV and E^ = Ep conjugated with 
the particle numbers Nr and Nl- This initial state de- 
termines the bare Keldysh Green functions. 

We will consider only the zero temperature limit. It is 
convenient to switch^ to the interaction representation 
Hq — ► Ho — E\Nr — E^Nl- This transformation induces 
a time dependence in the electron creation and annihila- 
tion operators. As the result each exponent in Eq. ([9]) is 
multiplied by exjp(ieVt[n^ + njj//i). 

In the Keldysh formulation^ the backscattering cur- 
rents (fTTj) . (fT2|) are evaluated as 

I b c % = (0|S(-oo, 0)j£S(0, -oo)|0>, (13) 

where |0) is the ground state for the Hamiltonian Hq, Eq. 
([7]), and S(t,t') the evolution operator for H' from t' to 
t in the interaction representation with respect to Hq. 
The result of this calculation depends on the elements of 
the matrix ([8]) , which describe the low-energy degrees of 
freedom and depend on the microscopic details. Several 
regimes are possible^ at different values of the parameters 
g s > 0, g c > 0, a and /?. In this paper we focus on one 
particular regime, in which the main contribution to the 
rectification current comes from backscattering operators 
17(1,0), U(0,-1) and U{-1,1). 

IV. RECTIFICATION CURRENTS 

In order to calculate the current (|13p we will expand 
the evolution operator in powers of U(n,m). Such per- 
turbative approach is valid only if U < Ep. The details 
of the perturbative calculation are discussed in Appendix 

o 

The currents (|13p can be estimated using a renormal- 
ization group procedure^. As we change the energy scale 



E, the backscattering amplitudes U{w[,n\) scale as 

U (n, to; E) - U(n,m)(E/E F ) z ^ m \ (14) 

where the scaling dimensions are 

z(n,m) = n 2 [g c (l + af + g s (l + (3) 2 ] +to 2 [ 5c (1 - a) 2 
+.9,(1 - I3) 2 } + 2nm[g c {\ - a 2 ) - g s (l - f3 2 )} - 1 

= n 2 A + m 2 B + 2nmC - 1. (15) 

The renormalization group (RG) stops at the scale of the 
order E ~ eV. At this scale the backscattering current 
can be represented as 7* s s = Vr c , s (V), where the effec- 
tive reflection coefficient r c s (V) is given by the sum of 
contributions of the formic 

(const )U(ni, mi; E = eV)U(ri2,m2; eV) . . . U (n p , to p ; eV) 

(16) 

Such a perturbative expansion can be used as long as 

U(n,m;E = eV) <E F (17) 

for every (n, m). This condition defines the RG cutoff 
voltage V* such that 

U{no, to ; E = eV*) = E F (18) 

for the most relevant operator U(riQ, mo). The RG proce- 
dure cannot be continued to lower energy scales E < V* . 

One expects that the leading contribution to the 
backscattering current emerges in the second order in 
U(n,m), if the above condition is satisfied. The lead- 
ing contribution to the rectification current may however 
emerge in the third order. Indeed, the second order con- 
tributions to the charge current were computed in Ref. 
[l4l . The spin current can be found in exactly the same 
way. The result is 

llf ] {V) ~ ^(const)|[/(n,TO)| 2 |T/| 2z (^ m ) +1 sign(F). 

If the (unrenormalized) U(n, m) were independent of the 
voltage, the above current would be an odd function of 
the bias and hence would not contribute to the rectifica- 
tion current. The backscattering amplitudes depend^ on 
the charge densities kp va though, which in turn depend 
on the voltage in our model^. The voltage-dependent cor- 
rections to the amplitudes are linear in the voltage at low 
bias eV -C Ep. Hence, the second order contributions to 
the rectification currents scale as U 2 \V\ 2z ^ n,m ^ +2 (see Ap- 
pendix [BJ . The additional factor of V makes the second 
order contribution smaller than the leading third order 
contribution (|23[) at sufficiently high impurity strength 
U *C Ep (as shown in Appendix [Bl U/Ep must exceed 
[V/Ep] 1+Z ^- Z ^-^ 1 -^). Note that the second or- 
der contribution to the rectification current is nonzero 
even for a symmetric potential U(x) and emerges solely 
due to the voltage dependence of the injected charge den- 
sity (cf. Sec. |TIJ). The leading third order contribution 
emerges solely due to the asymmetry of the scatterer. 
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The main third order contribution comes from the 
three backscattering operators, most relevant in the 
renormalization group sense (small z(n,m), Eq. (1151) ). 
They are identified in Appendix [Bj Under conditions 
(|B6IB7IB8IB12IB14IB15|) . the most relevant operator is 
[7(1, 0) , the second most relevant {7(0,— 1), and the 
third most relevant U(— 1, 1). The cutoff voltage V* 
is determined by the scaling dimension z(l,0), eV* ~ 
E F {U/E F fl^- A ^. The leading non-zero third order con- 
tributions to the spin and charge currents come from the 
product of the above three operators in the Keldysh per- 
turbation theory (see Appendix B). This leads to 



jbs ^ jj3y2(A+B-C-l) 



(20) 



This contribution dominates the spin rectification current 
at 



E F (U/E F ) 1/Vi+2CI ~ 2Bi = eV** >eV> eV* (21) 



as is clear from the comparison with the leading second 
order contribution I b s s 2 ~ U 2 V 2A , Eq. (fETO]) . Interest- 
ingly, the current ([20]) grows as the voltage decreases in 
the regime l|B6|B7|B8lB12|B14|B15jl . 



However, does the current ([20]) actually contribute to 
the rectification effect? In general, ([20[) is the sum of odd 
and even functions of the voltage and only the even part 
is important for us. One might naively expect that such a 
contribution has the same order of magnitude for the spin 
and charge currents. A direct calculation shows, however, 
that this is not the case and the spin rectification current 
is much greater than the charge rectification current. 

In order to calculate the prefactors in the right hand 
side of Eq. (j20|) one has to employ the Keldysh formal- 
ism. The details are explained in Appendix [C] Here 
let us shortly summarize the essential steps: The third 
order Keldysh contribution reduces to the integral of 
P{h,t 2 ,h) = (T c exp(# T (ti)+ieV^/l)exp(-%(i 2 ) - 
ieVt^/h) exp(i[-^(t 3 ) + <pi(t 3 )])) over (h - t 3 ) and 
{ts ~ £3)) where T c denotes time ordering along the 
Keldysh contour —00 — ► — * —00 and the angular brack- 
ets denote the average with respect to the ground state of 
the non- interacting Hamiltonian ((7|) . The integration can 
be performed analytically as discussed in Appendix [C] 



One finds 

rbs _ 16eT C 



irk 3 



sign(eV) 



eVr c 



a+b+c-2 



r(i-o)r(i-6) 



x r(2 - a - b - c)T(a + b - 1) sin — 

. irb . 7r(a + b) . , , 

x sm — sm sm 7r(a + b + c) 

x Re[{7(l, 0){7(-l, 1)17(0, -1)], 



(22) 



jbs 



16 77 na irb 
Trfi 2 2 2 



eVr c 



a+b+c-2 



xT(a + b- l)r(2 - a - b - c)T(l - a)T(l - b) 

7r(a + b + c) 



{lm[E/(l,0)Z7(-l,l)Z7(0,-l)] 



— 1 I cos ■ 

" 2 

, 7rc ir(a-b) 7r(a + b + c) irta + b) 

x sm 1- cos sm h sm 

1 2 2 2 2 

x cos ^( a +^ + c ) ] + lRe[[/(l, 0)U(-1, 1)17(0, -1)] 



7r(a — b) 
x sm _ 



sin7r(a + b + c)sign(eU) j, 



(23) 



where a = 2A - 2C*, b = 2B - 2C, c = 2C and r c ~ 
h/E F is the ultraviolet cutoff time. The charge current 
([2"2"| is an odd function of the voltage and hence does not 
contribute to the rectification effect. The spin current 
(|23p is a sum of an even and odd functions and hence 
determines the spin rectification current 



Tr 16-r 2 . ira . irb ir(a + b + c) 

L = — — sm — sm — cos 

s nh 2 2 2 2 



eVr c 



a+b+c-2 



xT(a+b- l)r(2 - a - b - c)T(l - a)T(l - b) 

x Im[C/(l, Q)C/(-l, 1)17(0, -l)][cos sin <^±±±A 



sm ■ 



. 7r(a + b) ir(a + b + c). 
sm cos 



(24) 



It is non-zero if Im[{7(l, 0)«7(-l, 1)U(0, -1)] ^ 0, which 
is satisfied for asymmetric potentials. The leading con- 
tribution to the charge rectification currents comes from 
other terms in the perturbation expansion. Thus, we ex- 
pect that in the region (|B6IB7IB8IB12IB14IB15|) . the spin 
rectification current exceeds the charge rectification cur- 
rent in an appropriate interval of voltages ([2"Tj) . The dif- 
ference between the spin and charge rectification current 
can be easily understood from the limit A = B. In that 
case the charge current changes its sign under the trans- 
formation {7(1,0) <-> {7(0,-1), V -> -V. Since {7(1,0) 
and {7(0,-1) enter the current only in the combination 
{7(1, 0){7(0, —1), this means that the charge current must 
be an odd function of the voltage bias. A similar argu- 
ment shows that at A = B the spin rectification current 
is an even function of the voltage in agreement with Eq. 
(US]). 

The voltage dependence of the spin rectification cur- 
rent is illustrated in Fig. O The expression ([2H)) de- 
scribes the current in the voltage interval V** > V > V* . 
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In this interval the current increases as the voltage de- 
creases in the regime (|B6IB7IB8IB12IB14IB15)) . At lower 
voltages the perturbation theory breaks down. The cur- 
rent must decrease as the voltage decreases below V* 
and eventually reach at V — 0. At higher voltages, 
Ep > eV > eV**, the second order rectification current 
(|19[) dominates. The leading second order contribution 
II ~ |J7(1, 0)| 2 y 2z ^ 1: °) +2 grows as the voltage increases. 
The charge rectification current has the same order of 
magnitude as the spin current. 

The Tomonaga-Luttinger model cannot be used for the 
highest voltage region Ep ~ eV. 

It is easier to detect charge currents than spin currents. 
However, the measurement of the spin current can be 
reduced to the measurement of charge currents: Let us 
split the right end of the wire into two branches and 
place them in opposite strong magnetic fields so that only 
electrons with one spin orientation can propagate in each 
branch. If both branches are grounded, they still inject 
exactly the same charge and spin currents into the wire as 
one unpolarized lead. However, the current generated in 
the wire will split between two branches into the currents 
of spin-up and spin-down electrons. If they are opposite 
then pure spin current is generated. 



APPENDIX A: HIGH POTENTIAL BARRIER 

In this appendix we first briefly consider the model 
of non-interacting electrons, Sec. [TTJ and then a simple 
Hartree-type model for strongly interacting electrons. 

1. Model without interaction 

We consider non-interacting electrons in the presence 
of the potential 

U[x) = uiS(x) +u%8{x - a). (Al) 

The transmission coefficient can be found from elemen- 
tary quantum mechanics, 



(1 — 2s\S2 sin 2 ka) 2 + (si + S2 + S1S2 sin2fca) 2 ' 

(A2) 

where E = H 2 k 2 /2m and Sj = mui/kh 2 . The spin and 
charge rectification currents can be computed from Eqs. 
(I1I2|) . Fig. [3] shows their voltage dependence for a certain 
choice of u%, U2, the voltage bias V and the magnetic field 
H. 



V. CONCLUSIONS 

In this paper, we have shown that rectification in quan- 
tum wires in a uniform magnetic field can lead to a spin 
current that largely exceeds the charge current. The pa- 
per focuses on the regime of low voltages and weak asym- 
metric potentials in which the perturbation theory pro- 
vides quantitatively exact predictions. Qualitatively the 
same behavior is expected up to eV, U ~ Ep . The spin 
rectification effect is solely due to the properties of the 
wire and does not require time-dependent magnetic fields 
or spin polarized injection as from magnetic electrodes. 
The currents are driven by the voltage source only. In 
an interval of low voltages the spin current grows as the 
voltage decreases. In contrast to some other situations, 
the z-component of the total spin conserves and hence 
the dc spin current is constant throughout the system. 
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2. Model with interaction 

It is difficult to find a general analytic expression for 
the current in the regime when both the electron interac- 
tion and potential barrier are strong. If all characteristic 
energies, U, eV, h 2 /[ma 2 ] and the typical potential en- 
ergy of an electron Ep, are of the order of Ep then one 
can estimate the spin and charge rectification currents 
with dimensional analysis: I r c ~ eEp/h, JJ ~ Ep. 

To obtain a qualitative picture of the interaction effects 
in the case of a high potential barrier (|A1[) . we restrict 
our discussion to a simple model in the spirit of the zero- 
mode approximation 2 ^. We assume that electrons move 
in a self-consistent Hartree-type field. In our ansatz the 
self-consistent field takes three different constant values 
Vl, Vm and Vr on the left of the potential barrier, be- 
tween two (5-function scatterers and on the right of the 
potential barrier. In the spirit of the Luttinger liquid 
model, we assume that the constants Vp, Vm and Vr 
are proportional to the average charge density in the re- 
spective regions, e.g., Vm = „ Jq dxp(x), where 7 is the 
interaction constant. 

A result is shown in Fig. [4j We see that the voltage 
dependence of the spin and charge rectification current 
exhibits a behavior similar to the non-interacting case. 

APPENDIX B: ESTIMATION OF HIGHER 
PERTURB ATIVE ORDERS 

In this appendix we compare contributions to the rec- 
tification currents from different orders of perturbation 
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theory. We focus on the regime when the third order 
contribution dominates. The appendix contains 5 sub- 
sections and has the following structure: 1) We introduce 
a parametrization for the scaling dimensions (fT5j) . 2) We 
discuss the operators most relevant in the RG sense. 3) 
We determine at what conditions the second order contri- 
bution to the rectification current dominates. Subsection 
3 also contains a lemma which is important in subsection 
4. 4) We determine at what conditions the third order 
contribution to the current dominates. 5) We estimate 
the voltages and currents at which the spin rectification 
current can exceed the charge rectification current in re- 
alistic systems. 

As shown in Refs. andd, there are two effects lead- 
ing to rectification in Luttinger liquids, which are here 
very shortly summarized: The density-driven and the 
asymmetry- driven rectification effects. The former ap- 
pears at second order in U. It appears because the 
backscattering potential depends on the particle densi- 
ties in the system, which in turn are modified by the 
external voltage bias. The leading order backscattering 
currents are of the formM I bs (V) ~ sign(F){7 2 |U| a so 
that the rectification currents, F = [I bs (V) + I bs {-V)} 
vanish. Due to the density dependence, however, an ex- 
pansion of U to linear order in V cancels the sign(V), 
and we obtain a rectification current V ~ [/ 2 |1/|" +1 . 

The asymmetry- driven rectification effect appears at 
third order in U. It is due solely to the spatial asym- 
metry of the potential U (x) : Due to backscattering off 
U, screening charges accumulate close to the impurity. 
Those create an electrostatic nonequilibrium backscat- 
tering potential W(x) for incident particles, leading to 
an effective potential U{x) = U(x) + W(x). The spatial 
distribution of charges follows from the shape of U(x) 
and the applied voltage bias. An asymmetric U(x) leads 
to different electrostatic potentials for positive or nega- 
tive bias, and hence to rectification. If we expand the 
current, I bs ~ U 2 ~ U 2 + UW + . . . , the asymmetry 
appears first at order UW. Since the charge density in 
the vicinity of the impurity is modified by the modifica- 
tion of the particle current through backscattering, W 
itself is (self-consistently) related to the backscattering 
current as W ~ I bs . Hence W ~ U 2 , so that the asym- 
metric rectification effect appears first at third order in 
U, I r ~ UW ~ U 3 . 

The main result of this paper are expressions for the 
currents that result from the perturbation theory at third 
order in the impurity potential U. In this appendix we 
show that the considered contribution indeed dominates 
the second and other third order expressions in the region 
defined by Eq. (|B6IB7IB8IB12IB14IB15|) . In addition, 
we give the proof that higher perturbative orders N > 4 
cannot exceed these values in the considered range of the 
system parameters g c ,g s ,a and (3. Unless we want to 
emphasize the correct dimensions, we set Ep = 1, e = 1 
and h = 1 in this appendix. We assume that U < Ep 
and eV < E F . 

An important observation is the following: In Eq. (TiTj| , 



Si n i — Yli m i = 0- This follows from the fact that in 
the absence of backscattering the numbers of right- and 
left-movers with different spin orientations conserve. 



1. Parametrization of scaling dimensions 



According to Eq. (JT| 

z(n, m) = n 2 A + m 2 B + 2nmC — 1, 



where 



A = [g c (l + a) 2 + g s (l + f3) 2 } 
B=[g c {l-a) 2 +g s (l-(3) 2 ] 
C=[g c (l-a 2 )-g s (l-f3 2 )} 



(Bl) 

(B2) 
(B3) 
(B4) 



Since g c and g a are positive, A and B are also positive. 
C can have any sign. It satisfies the inequality 



\C\< 



(B5) 



Indeed, AB - C 2 = ig c g s (l - a/3) 2 > 0. Any values of 
A,B>0 and -a/AB < C < \[AB are possible. For 
example, one can set a = (3 = (\/^4 — V^B) / (V^4 + vB), 
= (VA + VB) 2 [1 + C/VAB]/8, g s = (VA + Vb) 2 [1 - 




2. Most relevant operators 

Depending on the values of A, B and C many dif- 
ferent possibilities for relative importance of different 
backscattering operators U(m,n) exist. In the paper we 
focus on the situation when the most relevant operator 
is U(1,0), the second most relevant operator is U(Q, — 1) 
and the third most relevant operator is C7( — 1, 1) (cer- 
tainly, the scaling dimensions of the operators U(n,m) 
and U(—n, —m) are always the same). We will also as- 
sume that the operator C/(1,0) is relevant in the RG 
sense, i.e., 2(1,0) < 0. The analysis of the situation in 
which £7(0, — 1) is the most relevant operator, {7(1, 0) is 
the second most relevant and U(— 1, 1) is the third most 
relevant follows exactly the same lines. Similarly, little 
changes if {7(1, 1) is the third most relevant operator. 

The scaling dimensions of the three aforementioned op- 
erators are A—l, B—l and A+B - 2(7- 1. The following 
inequality must be satisfied in order for these operators 
to be most relevant backscattering operators: A — 1 < 
B — l < A + B — 2C — 1 < [all other scaling dimensions ] . 
Hence 



Since z(l,0) < 0, 



B > A > 2C. 



A<1. 



(B6) 



(B7) 



When are all other operators less relevant? We must 
consider three classes of operators: 1) {7(1, 1); 2) U(n, 0) 
and {7(0, n) with |n| > 1; 3) all other operators. 
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1) Since z(l, 1) = A + B + 2C - 1, one finds 

C > 0. (B8) 

2) z(0, re) = Bn 2 - 1 > z(n, 0) = An 2 - 1 > 4A - 1 > 
A + 5 - 2C - 1. Thus, 

3i|2C>B, (B9) 

3) z(n,m)-(A+.B-2C-l) = An 2 +BTO 2 + 2CreTO-(A+ 
B-2C)> An 2 + Bm 2 -C(n 2 + m 2 )-A-B + 2C = (A — 
C)(n 2 -1) + (B-C)(m 2 -1) > since B-C > A-C > 
in accordance with Eq. (jB6|) . |n|, |m| > 1 and either |n| 
or |rre| exceeds 1. Thus, case 3) gives no new restriction 
on A, B and C. 

3. Second order contribution to the current 

When is the second order contribution to the recti- 
fication current dominant? Any operator U{n, m) can 
be represented as [7(re, to) + VUi(n,m) + . . . , where 
Ui ~ U/Ef- Any second order contribution to the cur- 
rent which contains U only is an odd function of the volt- 
age bias. Indeed, any such contribution is proportional 
to U(n,rn)U*(n,m) = [7(re, to)[7(— re, — to). The trans- 
formation U «-» {/*, V — > — V changes the sign of the 
current. At the same time, the transformation U U* 
cannot change the second order current at all. Hence, it 
is odd in the voltage. The same argument applies to any 
perturbative contribution which contains only U, if every 
operator U(n,m) enters in the same power as U*(n, m). 
In particular, if only two operators U(n, 0) and £7(0, m) 
and their conjugate enter then the resulting current con- 
tribution is odd. 

Thus, all second order contributions to the rectification 
current must contain U\. As is clear from Eq. (fT9|) , the 
leading second order contribution is proportional to the 
square of the most relevant operator, \U(1, 0)| 2 . It scales 
as 

h ~ VU 2 V 2z ^ +1 ~ U 2 V 2A . (BIO) 

In this subsection we discuss at what conditions this 
contribution dominates for all V > V* (see Eq. ([15])), 
Since {7(1, 0) is the most relevant operator, its renormal- 
ized amplitude U(1,0;E = V) exceeds the renormalizcd 
amplitude of all other operators on every energy scale. 
At the same time it remains lower than 1 (i.e. Ep) for 
V > V*. This certainly means that the renormalizcd 
amplitudes are smaller than 1 for all other operators too. 
Hence, the product of any operators is smaller than the 
product of any two of them and that product cannot ex- 
ceed U 2 (1,0; E). This guarantees that the second order 
current (|B10[) exceeds any second or higher order contri- 
bution which contains any operator VUi(n, to). Thus, we 
have to compare I2 with higher order contributions to the 
rectification current which contain U only. Every such 
contribution is at least third order and contains at least 



one operator less relevant than £7(0,— 1) [if it contains 
f7(±l,0) and {7(0, ±1) only then it must contain VU\ as 
discussed above]. Thus, any rectification current contri- 
bution with U only cannot exceed u^y 2z { 1 fi)+ z ( 1 - 1 )+ 1 _ 
Comparison with (|B10|) at V ~ V* leads to the condition 

B>2C+1. (Bll) 

4. Third order contribution to the current 

The most interesting question is different. When 
does the third order contribution dominate the rectifi- 
cation current? We will focus on the third order con- 
tribution I 3 proportional to {7(1, 0){7(0, — 1){7(— 1, 1) at 
V ~ V* . Note that this contribution is proportional to 
~ V 2 (^+- B_C ' _1 J and hence scales as a negative power of 
the voltage, if 

A + S<C + 1. (B12) 

At V~ V*, U^V 1 ~ A . Thus, 

I 3 (V = V*)~V 2B - A - 2C+1 . (B13) 

We need to compare J3, Eq. (|B13[) . with the follow- 
ing types of contributions: 1) those containing at least 
three different operators [we treat a pair of U (n, m) and 
U(—n,—m) = U*(n,m) as one operator]; 2) those con- 
taining only one type of operators; 3) those containing 
two types of operators. 
Cases 1) and 2) are easy. 

1) Is contains the product of the three most rele- 
vant operators and hence always exceeds the product of 
any other three different operators at any energy scale 
Ep > V > V* . Any contribution with three different op- 
erators is the product of three different operators times 
perhaps some other combination of operators which can- 
not exceed 1 at Ep > V > V* . Hence it is smaller than 
h- 

2) Any contribution to the rectification current with 
only one type of operators must contain VU\. As dis- 
cussed in the previous subsection, the leading contribu- 
tion of such type emerges in the second order. It is I2, Eq. 
(iBTOlt . At V - V*, h(V = V*) ~ V 2 . The condition 
h (V*) < h(V*) means that 

2B < A + 2C +1. (B14) 

3) We have to consider three possibilities: 3.1) one 
operator has the form U (n, 0) and the second operator 
has the form U(k,m), to 7^ or one operator has the 
form [7(0, to) and the other one has the form {7(n, k), 
n ^ 0; 3.2) both operators have the form U(ni, 0) or both 
operators have the form [7(0, m*); 3.3) both operators 
have the form U(rii, TOj) with rii, TOj =/= 0. 

3.1). Let us assume that one operator has the form 
[7 (re, 0) and the second one is U(k,m). The case of the 
operators [7(0, to) and U(n, k) can be considered in ex- 
actly the same way. We must have the same number 
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of operators U(k,m) and U(—k,—m) in the perturba- 
tive contribution since the sum of the second indexes ±m 
must be 0. [The other cases are covered in 3.3).] From 
the analysis of the sum of the first indexes one concludes 
that the operators U (n, 0) and U(—n, 0) also enter in the 
same power. It follows from the previous subsection that 
the perturbative contribution must contain at least one 
Ui operator and hence is smaller than l^. Hence, it is 
also smaller than J3. 

3.2) We will focus on the case when both operators 
have the form U{ni,0). The case when both operators 
have the form V{0, rrii) is very similar and does not lead 
to a new restriction on A, B and C . The scaling dimen- 
sions of the operators U(n,0) are An 2 — 1. Operators 
with greater n are less relevant. Since the contribution 
contains two different operators, it must be at least third 
order [we treat U(n, 0) and U(— n, 0) as the same opera- 
tor!]. At least one of the two operators must have \rii\ > 1 
(otherwise all operators are t/(±l,0)). Thus, the contri- 
bution cannot exceed U 2 {l,0;E = V)U(2,0;E = V) ~ 
U 3 V 6A - 2 . The comparison with I 3 ~ ij3 V 2A+2B-2C-2 
at E F > V > V*, yields: 

B<2A + C. (B15) 

Note that the above condition is stronger than (|B9[1 . 

3.3) This case is easy: the contribution must be at 
least third order again. Both operators U (rii, mi) are less 
relevant than [7(1,0) and U(0, —1) and no more relevant 
than {7(1, —1). Thus, the contribution is automatically 
smaller than ^3 at any energy scale Ep > V > V* . 

We now have a full set of conditions at which the third 
order contribution dominates at V ~ V* and the spin 
rectification current scales as a negative power of the 
voltage. These are equations (|B6IB7IB8IB12IB14IB15|) . 

The above analysis shows that ^3 exceeds any contri- 
bution to the spin rectification current which does not 
contain VU\ in the whole region Ep > V > V* . 1% 
dominates the remaining contributions for any V > V*. 
The contributions become equal, I2 = I3, at V = V** — 
j/i/[2+2C-2B] < In thc i nterva i f voltages V** > V > V* , 
the spin rectification current is dominated by I3. At 
V > V**, the spin and charge rectification currents are 
dominated by 



5. Numerical estimates 

In order to get a feeling about the magnitude of the 
effect, let us consider a particular choice of parameters 
A = B = 7/12, C = 7/24, eV ~ 0ME F , eV* ~ 10~ 4 i; F . 
For such A, B and C the scaling dimensions of the three 
most relevant operators are the same. The inequali- 
ties (|B7IB8IB12IB14IB15P are satisfied. The equality 
A = B = 2C corresponds to a limiting case of (|B6[) . One 
finds that U - OMEp and eV** ~ 0AE F . Repeating 
the arguments of the previous section one can estimate 
the leading correction to -[3 as SI ~ (eV/Ep) 7 / 12 !^ <C -[3. 



The spin rectification current is the difference of two op- 
posite electric currents of the spin-up and -down electrons 
times H/[2e]. Even if E F is as low as ~ 0.1 meV, this still 
corresponds to the voltage V of the order of microvolts 
and thc currents^ ([24]) of spin-up and -down electrons of 
the order of picoamperes, i.e. within the ranges probed in 
experiments with semiconductor heterostructures. Cer- 
tainly, the current increases, if E F or V* is increased. 

APPENDIX C: EXPLICIT EVALUATION OF 
THE THIRD ORDER CURRENTS 

The charge or spin currents in the third order in the 
potentials U are evaluated from thc following perturba- 
tive expression: 

x (r„t;(rt t ,nj;0)&(m,,i7.j;« 1 )Er(I,,l i; ( 2 )), (CI) 

where the sum runs over indices satisfying n a +ma-+la- = 
for a =|, 4, Ck is the Keldysh contour —00 — > — > 00, 
T c the time order on Ck, and we omitted a constant 
prefactor. The operators U are given by 

U{n hni ;t) = U{n h n l )e t( - n ^ +n ^ teV/h e m ^^ t)+m ^i {t) . 

(C2) 

The most relevant expressions are those arising 
from thc combinations U(l, 0)U(0, — l)U(— 1, 1) and 
[/(-1,0)[/(0,1)[/(1,-1) (see Appendix[B]). The third or- 
der contributions to the current contain correlation func- 
tions of the form 

P(ti,t 2 ,t 3 ) = (T c e ±i[ ^(*i)-^(*=)-^T(t3)+<Ai(f3)]^ ( C3 ) 

We evaluate the correlation functions within thc 
quadratic model described by Eq. (JT]) and use the rela- 
tions © and (4> c , s (ti)4> c ,s(ti)) =-21n(i(ti-t 2 )/T c + a), 
with an infinitesimal S > and t c ~ h/E F the ultraviolet 
cutoff time. This leads to 

P(h,t 2: t 3 ) = (iTcih-t^/Tc + s) 20 - 2 * 

x (»T c (i a - h)/r c + S) 2 °~ 2B \iT c {U - t 2 )/r c + 6)~ 2C , 

(C4) 

where T c (tj— tj) — (ti—tj), if time stays later than tj on 
the Keldysh contour, and otherwise T c (tj — tj) — (tj — ti). 
The expression (|CJ4|) is independent of the ± signs in 
Eq. (IC31 . 

The spin and charge current contributions, propor- 
tional to 17(1, 0)1/(0, — 1)U(— 1, 1), arc complex conju- 
gate to those proportional to U(— 1, 0)U(0, l)U(l, —1) = 
U*(1,0)U*(0,-1)U*(-1,1). Thus, it is sufficient to cal- 
culate only the contributions of the first type. In the case 
of the charge current, their calculation reduces to the cal- 
culation of the following two integrals over the Keldysh 
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contour: 

J dtidt 3 P(*i,Q 1 t3)exp(ieV r ti/n) (C5) 

and 

J dt 2 dt 3 P(0,t 2 ,t 3 ) exp(-ieVt 2 /h). (C6) 

One of the times t\ and t 2 is zero since the current opera- 
tor is taken at t = in Eq. ([TU]) . The two integrals can be 
evaluated in exactly the same way. We will consider only 
the first integral. We find 8 integration regions. They 
correspond to 2 x 2 — 4 possibilities for the branches of 
the Keldysh contour on which t\ and £3 are located and 
two possible relations |ti| > | ^3 j or \t 3 \ > In all 8 
cases, we first integrate over t 3 . The integral reduces to 
the Euler B- function. Then we integrate over t±. This 
yields a T-function. Finally, we obtain Eq. (12"2")) . 



The spin current contains three contributions propor- 
tional to U(l,0)U(Q,-i)U(-l, 1). Two of them reduce 
to the integrals (|C5j) and (|C6j) . The third contribution is 
proportional to 



J dtidt 2 P(h,t 2 ,0) exp(ieV[t! - t 2 ]/h). (C7) 



Again we have eight integration regions determined by 
the choice of the branches of the Keldysh contour and 
the relations > \t 2 \ and \t 2 \ > \ti\. In each region 
it is convenient to introduce new integration variables: 
t = \t\ — t 2 \ and t = miii(ti,t 2 ). The integration over t 
reduces to a i?-function. The integration over r produces 
an additional T- function factor. Finally, one obtains Eq. 

(EH). 
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FIG. 1: Sketch of the one-dimensional conductor connected 
to two electrodes on both ends. Currents are driven through 
a voltage bias V that is applied on the left electrode while the 
right electrode is kept on ground. The system is magnetized 
by the field H. Electrons are backscattered off the asymmetric 
potential U(x). U(x) 7^ in the region of size au ~ 1/kF- 
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FIG. 2: Double-well potential with quasistationary levels. 
The transmission coefficient is maximal in the shaded regions. 
The narrow potentials u\(x) and u%{x) are centered at the 
positions x — and x — a (a < kp 1 ), respectively, and are 
modeled by 8- functions in Eq. (|Aip . 



13 



0.15 - 



0.00 



-0.15 - 



-0.30 - 




0.0 0.2 0.4 0.6 

V/V 



FIG. 3: Normalized charge rectification current I c /Ic.o and 
spin rectification current I a /I 3 o versus applied voltage V/Vo 
for non-interacting electrons with Ef — 400eo, \xH = 75eo, 
ui — 50eocs and ui = — 50eoa, where eo = h 2 /ma 2 (see Fig. [2] 
and Appendix|Aj • IcO = 50eeo /fr, IsO = 25eo and Vo = 50eo / e 
are arbitrary reference currents and voltage. 
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FIG. 4: Normalized charge rectification current I c /Ic0 and 
spin rectification current I a /I 3 o versus applied voltage V/Vo 
for interacting electrons with Ef = lOOeo, liH = 25eo, 7 = 
12.6eoa/e, ui = 25eoa and U2 = 50eoa, where eo = fi 2 /ma 2 
(see Fig. |2] and Appendix |A|)- IcO = 50eeo/ft, ho = 25eo, and 
Vo — 50eo/e are arbitrary reference currents and voltage. 
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FIG. 5: Qualitative representation of the spin rectification 
current. The spin current exceeds the charge current and 
follows a power-law dependence on the voltage with a negative 
exponent in the interval of voltages V < V < V** . 



